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Abstract 

We study three convolutions of polynomials that are inspired by free probability. We 
define these to be the expected characteristic polynomials of certain random matrices. The 
symmetric additive and multiplicative convolutions have been studied for a century. The 
asymmetric additive convolution, and the connection of all of them with random matrices, 
appears new. We prove that these convolutions produce real rooted polynomials and provide 
strong bounds on the locations of the roots of these polynomials. 
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1 Introduction 

We study three convolutions on polynomials that are inspired by free probability. Instead 
of capturing the distributions of eigenvalues of random matrices, these capture the expected 
characteristic polynomials of random matrices. Two of the convolutions have been clasically 
studied. The third, and the connection of all of them with random matrices, appears new. 

For a matrix M, we define 

X:r (M) = det(x/- M) 

to be the characteristic polynomial of the matrix M in the variable x. 

The convolutions we study may be defined in terms of integrals over the group of orthonormal 
matrices. For d-dimensional symmetric matrices A and B with characteristic polynomials p and 
q, we define the symmetric additive convolution of p and q to be 

p{x) Eld q{x) =Exx{A + QBQ^) , (1) 

where the expectation is taken over orthonormal matrices Q sampled according to the Haar 
measure. 

When A and B are positive semidefinite, we define we dehne the symmetric multiplicative 
convolution of p and q to be 

p{x) [x]^ q{x) = E Xx {AQBQ^) , (2) 

let' 

where the expectation is taken over random orthonormal matrices Q sampled according to the 
Haar measure. 

For arbitrary d-dimensional square matrices A and B that are not necessarily symmetric 
such that AA^ and BB^ have characteristic polynomials p and q, we dehne the asymmetric 
multiplicative convolution of p and q to be 

p{x) EEd q{x) = E Xx ((^ + RBQ){A + RBQ)'^) , 

R,Q 

where the expectation is taken over random orthonormal matrices R and Q sampled according 
to the Haar measure. 

We derive expressions for these expected characteristic polynomials, prove that all of the 
resulting polynomials are real rooted, and prove bounds on the locations of their roots. 

Free probability theory (see [Voi97t INSDBt lAGZlO] considers the limits for large matrices of 
the moments eigenvalues of matrices such as those above. For each operation, free probability 
provides a transform of the moments of the eigenvalues of the individual matrices from which 
one can easily derive the transform of the moments of the eigenvalues of the resulting matrices. 
We show that for our “finite free convolutions” these transforms provide upper bounds on the 
roots of the resulting polynomials. For the symmetric additive convolution, this transform is 
equivalent to the barrier function studied in |BSS12l IMSSlSb] . In |MSS15cj , we use the result 
of the present paper to prove the existence of bipartite Ramanujan graphs of every number of 
vertices. 
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1.1 Summary of Results 

We derive the following expression for the symmetric additive convolution of two polynomials. 

Theorem 1.1. Let 


= *(—l)*aj, and q{x) = '^^x'^ *(— 


i=0 


i=0 


Then, 


p{x)mdq{x) = Y,x’^ ^(-1)^ Z] 


k=0 


i-\-j=k 


d\{d -i-j) 


This is equivalent to the identities 


■id 1 ^ 

p{x)mdq{x) = — ^(d-i)!6ipW(x) = -j^^p^"\x)q^'^~^\y) 


i=0 


dl 


i=0 


(3) 


where denotes the ith derivative of p. 

The convolution ([3]) was studied by Walsh |Wal22| , who proved results including the following 
theorem. 


Theorem 1.2. If p and q are real rooted polynomials of degree d, then so is q. Moreover, 

maxroot q) < maxroot (p) + maxroot {q) . 

See also |Mar66| and |RS02l Theorem 5.3.1]. 

In Theorem 11.71 we strengthen this bound on the maximum root. Our result is much tighter 
in the case that most of the roots of p and q are far from their maximum roots. 

We derive the following expression for the asymmetric additive convolution. 

Theorem 1.3. For 

d d 

p{x) = ^ q{x) = ^ l)*6j, 

having only nonnegative real roots, 

k^—Q z~\~j^—k 

This is equivalent to the identity 

/ 1 \ ^ ^ 

p{x) ESrf q{x) = ^((d - iy?)\{DxDfp{x). 

We are not aware of previous studies of this convolution. We prove the following theorem in 
Section [3l 
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Theorem 1.4. If p and q have only nonnegative real roots, then the same is true of pSEEB^ q. 

Theorem 1.5. For 


d d 

p{x) = {—ly Qj and q{x) = '^^x'^~'‘{—lYbi 

having only nonnegative real roots, 

d , 

p{x) ma qix) = x'^-^i-iy^. 

i=0 vj 

This convolution was studied by Szego |Sze22j . who proved the following theorem. 

Theorem 1.6. If p and q have only nonnegative real roots, then the same is true of pEKld q. 
Moreover, 

maxroot {plEld q) < maxroot {p) maxroot {q). 

We strengthen this result in Theorem 11.81 


1.2 Transforms 

In free probability, each of these three convolutions comes equipped with a natural transform of 
probability measures. We define analogous transforms on polynomials and use them to bound 
the extreme roots of the convolutions of polynomials. 

We will view a vector (Ai,..., A^) as providing a discrete distribution that takes each value 
Aj with probability \/d. The Cauchy/Hilbert/Stieltjes transform of such a distribution is the 
function 

1 1 
2=1 

Given a polynomial p with roots Ai,..., A^, we similarly define 

Qp (x) = G\ (x). 

We will prove theorems about the inverse of the Cauchy transform, which we define by 

KEp {w) = max{x : Gp {x) = tc} . 

For a real rooted polynomial p, and thus for real Ai,..., A^, it is the value of x that is larger 
than all the Aj for which Gp {x) = w. 

As Gp {x) = 3 ^^, 


Gp (x) = w 


P(x) - iP'(x) = 0 . 

wd 


This tells us that 


KEp (w) = maxroot (C/i/^dp) , 
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where we define 


Uap{x) = p{x) — aDp{x). 

Voiculescu’s R-transform of the probability distribution that is uniform on A is given by 
TZx (w) = JC\ (w) — 1/w. We overload this notation to produce a transform on polynomials 

TZp [w) = Kp (w) — Ijw. 

If A and p, are probability distributions on the reals, and if A ffl ^ is their free additive 
convolution, then 

(w) = 7^A (w) + Up, (w) . 

For our finite additive convolution, we obtain an analogous inequality whose proof may be found 
in Section O 

Theorem 1.7. For degree d polynomials p and q that have only real roots, 

(^) ^ '^p ('^) + (^) ■ (^) 

We will often write this inequality as 

maxroot (Uaip ffld ?)) + da < maxroot (Uap) + maxroot {Uaq ). 

To define the transform for the asymmetric additive convolution, we define S to be the oper¬ 
ation that converts a polynomial p{x) into the polynomial p(x^). If p has only positive real roots 
Aj, then §p has roots ±x/Xi. We overload this notation so that if A is a probability distribution 
supported on the nonnegative reals, then §A is the corresponding probability distribution on 
±v^. 

If A and p are probability distributions on positive numbers, and if A I++I p is their free 
asymmetric convolution, then 

^SAEUSp (^) = (w) + (w) . 

We obtain an analogous inequality whose proof appears in Section 14.31 
Theorem 1.8. For degree d polynomials p and q having only nonnegative real roots, 

(^) ^ (^) + (w) . 

To bound the roots of the finite multiplicative convolution, we employ a variant of Voiculescu’s 
S-transform. We first define a variant of the moment transform, which we write as a power series 
in l/z instead of in z: 

Mpiz) = zgp{z)-l = ^Y,Y.(^'\ . 

i=l j>l ^ ^ 

For a polynomial p having only nonnegative real roots and a ^ > 0, 

z > maxroot (p) M.p {z) < oo. 

We define the inverse of this transform, M.p {w), to be the largest 2 so that M.p (z) = w, and 

This is the reciprocal of the usual S'-transform. We prove the following bound on this 
transformation in Section [4.21 
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Theorem 1.9. For p and q having only nonnegative real roots and aw>0, 


(^) ^ (^) (^) • 


1.3 Notation 

Let ¥{d) be the family of real rooted polynomials of degree d with positive leading coefficient, 
and let P be the union over d of P((i). Let P^((i) be the subset of these polynomials having only 
nonnegative roots. We let P"*" be the union of P^((i) over all d > 1. We also define P“((i) and 
P“ to be the set of polynomials having only nonpositive roots. 

For a function f{x), we write the derivative as Df{x). For a number a, we let Ua be the 
operator that maps / to / — aDf. That is, Ua is multiplication by 1 — aD. 

2 Convolution Formulas 

In this section, we will derive explicit formulas for the additive (both symmetric and non- 
symmetric) and multiplicative convolutions of two polynomials in terms of their coefficients. 

We denote the coefficient of (—of the characteric polynomial of a d-dimensional 
matrix A by ek{A), which we recall is the kth elementary symmetric function of the eigenvalues 
of A. For a subsets S', T C {1,..., d}, we write T(S, T) for the submatrix of A indexed by rows 
in S and columns in T. We write ^(S,:) for the submatrix containing rows in S. We recall that 

ei(A) = det{A{S, S)). 

\S\=k 


2.1 Multiplicative Convolution 

We establish Theorem 11.51 bv establishing the identity for each of the coefficients individually. 


Theorem 2.1. For d x d symmetric matrices A and B and random orthonormal Q: 


Eae,(AQBQ^) = 

\k) 


Proof. We first show that we can assume without loss of generality that A and B are diagonal. 
To see this, diagonalize these matrices as A = UCU'^ and B = VDV"^ where U and V are 
orthonormal. We then have that 


Eg det(x/ - AQBQ^) = Eg det(x/ - UCU^QVDV'^Q'^) 

= Eg det(x/ - CU^QVDV^Q^U) 
= Eg det{xl -CQDQ^), 

as U^QV is distributed according to the Haar mesure. 
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Assuming that A and B are diagonal, we compute 


EgekiAQBQ'^) =E ^ det{A{S,:)QBQ{S,:f) 

\S\=k 

= Eq det(A(5, S)Q{S, :)BQ{S, :)'^) since A is diagonal 

\S\=k 

= EqY^ det{A{S,S))det{Q{S,:)BQ{S,:f) 

\s\=k 

= EqEp ^ det{A{S,S))det{{PQ){S,:)B{PQ){S,:f) 

\s\=k 

where P is a random permutation matrix, since PQ '=* Q 


EqEp det{A{S,S))det{P{S,:)QBQ^P{S,:y^ 
\s\=k 


= Eq ^ det(A(S,S)) (Epdet(P(P,:)QPQ^P(P,:f)) 
\\s\=k j 

since the distribution of P{S ,:) is independent of S 
= EQefc(A) ( ^ ^ det((QPQ^)(T,r) 

\ Uj \T\=k 


= E, 


Q- 


ek{A)ekiQBQ^ 


it) 

ek{A)ek{B) 

it) ■ 


□ 


2.2 Symmetric Additive Convolution 

We now prove Theorem II.II bv establishing it for each coefficient. 

Theorem 2.2. If A and B are d x d symmetric matrices then 

EQBk{A + QBQ'^)= Y 

i+j=k ' 

where Q is a random orthogonormal matrix. 

Proof. Writing A = UCU'^ and B = VDV^ where C and D are diagonal and U and V are 
orthonormal, we derive 

EQefc(A + QBQ'^) = EqekiC + U^QVDV^Q'^U) = Eq^kiC + QDQ^), 
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since U"^QV is a random orthonormal matrix. Thus, we may assume that A and B are diagonal. 
Given this assumption, we let e* denote the elementary unit vector in direction i and compute 


EgekiA + QBQ^) = Egcfc ^ aieiej + ^ hi{Qei){Qeif j 

i<d j 

= Eq ^ ek i'^aieicf+ '^bi{Qei){Qeif\ 
5,T:|5| + |r|=fc \ies ieT / 

= Eq ^ efc {Mst{A{S, S) © B{T, r))Mj^) 

\S\ + \T\=k 

= Eq ^ det{M^TMsTiA{S, S) © B{T, T)) 

\S\ + \T\=k 

= clet A{S, S) det B{T, T) ■ (Eg det(MjrMsr)) 

|S| + |T|=fc 


where Mst = 'S')| ((5/)(:, T)] is d x k with columns {ei}i^s and {Qei}i^T- 

It now turns out that the expectation Eg det(Mj,pM 5 'r) depends only on |5| =: i 
|r| =: j. Fixing S and T and removing subscripts to ease notation, we have 


EQdet(M^M) = EQe,(/(:,5)/(:,5f )e,(/(:,5f g/(: r)/(:, T)^g^/(:, 5)) 

= EQEn(l) • e,(/(:,5)'^ng/(: r)/(:, r)^g'^n^/(:, 5)) 

where 11 is a random permutation matrix, since Iig Q 
= EgE^e,•(/(:, RfQ{:,T)Q{:, Tfl{R ,:)) 

where i? is a uniformly random subset of [d] of size \S\ = d — i 

= ^QJ^ E writing P = g(:,T)g(:,rf 

\d-i) \R\=d-i 

“i^OT^-r E E 

\d-i) \R\=d-iWcR,\W\=j 


= E, 


Q / d 


E 


/_g^.)e,(P(,r.w.)) 


(d-i) Wcld],\W\=j 
since each W of size j = k — i<d — i appears 

= EQA^e,(P) 


d-j 

d-i-j 


times 


idt.) 




(/J 

(d-i)!(d-j)! 

d\{d-k)\ 


■ 1 since P is a projection of rank j. 


and 



Thus we may write the quantity of interest as 

S E det(A(S,S))det(B(T,T)) 

i+j=k '' \S\=i,\T\=j 

as desired. 


E 

i-\-j=k 


{d-i)\{d-j)\ 

d\{d-k)\ 


ei{A)ej{B), 


□ 


We remark that the above formula may also be obtained using the methods (based on estab¬ 
lishing rotational invariance) we use in the next section for the asymmetric additive convolution. 


2.3 Asymmetric Additive Convolution 

Recall that the asymmetric additive convolution of the characteristic polynomials of matrices A 
and B is defined as: 

h[A, B]{x) = EQ^ijdet(x/ — (A -|- QBR^){A + QBR’")'^) = p{x) 1 ++ !^ q{x), 

where the expectation is taken over random orthonormal matrices and 

P{x) = Xx , q{x) = Xx (BB^) . 

In this section, we will derive the explicit formula: 

h[A,B]{x) = ^x‘^->^{-lf (d^-k)l ^ T 

Observe that both sides of the equation are invariant under left and right multiplication of A 
and B by orthonormal matrices, so by taking singular value decompositions it suffices to prove 
it for the case when A and B are diagonal. This will be accomplished in two steps. First, we will 
pass from random orthonormal matrices to random signed permutations. A signed permutation 
matrix is a matrix whose nonzero entries are ±1 that has exactly one nonzero entry in each row 
and column. A uniform random signed permutation matrix may be obtained by multiplying a 
random diagonal with uniform ±1 entries by a uniform random permutation matrix. 

Define 

f[A,B]{x) Ep^sdet{xl - {A + PBS^){A + PBS^f), 
where P and S are uniform random signed permutation matrices. 

Theorem 2.3 (Asymetric Additive Quadrature). For all d dimensional square matrices A and 

B, 

f[A,B]ix) = h[A,B]ix). (7) 

This theorem establishes a quadrature result: it shows that the integral over the the or¬ 
thonormal group may be computed by taking the expectation over a discrete subset of that 
group. Using similar techniques, we may establish the following theorem. 

Theorem 2.4 (Symmetric Additive Quadrature). For all d dimensional symmetric square ma¬ 
trices A and B, 

^Xx{A + QBQ^) =ExxiA + PBP^) , (8) 

where the first expectation is over orthonormal matrices and the second is over signed permuta¬ 
tion matrices. 
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Then, we will explicitly calculate the polynomial f[A,B](x) in the diagonal case. 


Theorem 2.5. When A and B are diagonal, 

nA,B](x] = ^ 

k —0 —k 

The combination of Theorems 12.31 and 12.51 immediately establishes Q . 

2.3.1 Quadrature 

Theorem 12.31 follows easily from the following lemma, which says that / is invariant under left 

and right multiplication oi A, B by orthonormal matrices. 

Lemma 2.6. For any A,B and orthonormal Q,R: 

f[A,QBR'^] = f[A,B] 

Proof of Theorem 1 3., '91 

h[A, B] = Eq^rxHA + QBR^){A + QBR^ f) 

= + PQBR^S'^){A + PQBR'^S^f) for any fixed P, S 

= Eq,rEp,5x((^ + PQBR^S^){A + PQBR^S^f) 

= EQ^Rf[A,QBR^] 

= ^Q,Rf[A, B] by Lemma [2T] 

= f[A,B]. 


□ 

We prove Lemma 12.61 by showing that it holds for orthonormal matrices acting on two- 
dimensional coordinate subspaces, and then appealing to the fact that such matrices generate 
all orthonormal matrices. For s ^ t, let Rst{Q) denote the d x d matrix whose nonzero entries 
are given by 

1. = Rsti9)[t,t] = cos(d) 

2. Rst(9)ls,t] = -Rst(9)lt,s] = sm(9) 

3. Rst(9)lk, fe] = 1 for A: / s, t 

The key observation is that any expected characteristic polynomial over a distribution that is 
invariant under a small subset of such rotations must be invariant under all of them. 

Lemma 2.7. Let A,B be d x d matrices and let s,t £ [d] be two distinct indices. Let P he a 
random matrix supported on 0{d) such that the distribution of P is the same as the distribution 
of PRstijT^/‘2,) for j = 0,1, 2,3. Then, for every angle 9, 

Ep det(A + PBP'^) = Ep det(A + PRst{9)BRst{9f P^). 
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Proof. The hypothesis tells us that 

3 

4 • Ep det(^ + PRst{e)BRst{efP^) = Ep ^ det{A + PRst{j7r/2)Rst{B)BRst{9f Rst{j7r/2f P^) 

j=0 

3 

= Ep ^ det(A + PRstijTT/2 + e)BRst{j7r/2 + OfP^) 
j=0 
3 

= Ep ^ det(i2,t(i7r/2 + B f P^APRst{j7r/2 + B) + B) 

j=0 
3 2 

= Ep EE (.Pgki{j'iTl2+e) 

j=0k=—2 

by Lemma 12.81 for some depending on P 

= Ep ^ cf + i^ + 

k=-2 

= Ep«), 

since the inner sum vanishes for \k\ = 1,2. 

Since this quantity is independent of 0, we may take 0 = 0 and the claim is proved. □ 


Lemma 2.8. Let A,B be d x d matrices and let 

det {Rst{B)ARst{Bf + B)=J2 Cfce'"® 

k 


Then = 0 for \k\ > 3. 


Proof. Recall that all 2 x 2 rotations may be diagonalized as 


Re = 


cos B 

— sin0 

= U 


0 ■ 

sin0 

cos B 

0 





where 


u = -^ 
V2 


1 

—i 


is independent of B. This implies that Rst{B) = VDV^ for diagonal D containing and e 
in positions s and t and ones elsewhere, with V independent of B. Thus, we see that 


det {RstiB)ARst{Bf + B) = det {RstiB)A + BRM) 

= det (VDV^A + BVDV^^ 

= det (^DV^AV + V^BVD^ 

Notice that the matrix M = DV^AV + V^BVD depends linearly on and that the 

(resp. e“*®) terms appear only in the sth (resp. tth) row and column of M. Since each monomial 
in the expansion of the determinant contains at most one entry from each row and each column 
and e*® • = 1, this implies that no terms of degree higher than two in e*® or appear. □ 
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We are now in a position to complete the proof of Lemma 12.61 and thereby Theorem 12.31 

Proof of Lemma \2.b\ Since each Q,i? is a product of a finite number of orthonormal matrices 
acting on two-dimensional coordinate subspaces, it suffices to show that 

f[A,Rst{e)DstB\ = f[A,B] = f[A,BDlRst{ef], 

for every rotation Rst{0) and diagonal reflection Dgt with D{s,s) = ±l,D(t,t) = ±1, and 
identity elsewhere. 

We will show the first equality. Since the distribution of P is invariant under reflections, it 
is immediate that 

f[A,DstB] = f[A,B], 

so replacing B by DgiB it suffices to show that 


f[A,Rst{e)B] = f[A,B] 


for every 6. For any (not necessarily symmetric) dxd matrix A, let dil(A) denote the symmetric 
2d X 2d matrix 


dil(A) 


■ 0 A 
A^ 0 ’ 


and let ^4 0 i? denote the 2d x 2d block matrix 


A 0 
0 B 


Observe that 


f[A,B]{x‘^) =Ep^sdetixI - dil{A) - {P (B S)d[\{B){P (B Sf). 


whence 


f[A, RMB]ix^) = Ep,sdet(x/ - dil(^) - (P 0 S)dil(P,t(0)P)(P 0 Sf) 

= Ep,s det(x/ - dil(A) - (P 0 S){Rst{0) 0 /)dil(P)(P,i(0) 0 ifP 0 sf) 

For every fixed S, the random matrix P0S' is invariant under right multiplication by Rst{j'^ /2f)(B 
I for j = 0,1, 2,3, since these are themselves signed permutation matrices. Thus, we may apply 
Lemma [27n to conclude that for every S', 


Epdet(x/ - dil(^) - (P 0 S){Rst{e) 0 /)dil(P)(P,t(0) 0 /)^P 0 S)^) 
= Ep det(x/ - dil(kl) - (P 0 S)dil(P)P 0 S)^). 


Taking expectations of both sides with respect to S yields f\A,Rst{0)B] = f[A,B], as desired. 
The argument for the other side is identical, except we consider KBRstf) instead of Rst{d)(B 

I. □ 
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2.3.2 Calculation for Diagonal Matrices 

Proof of Theorem \2 . 51 Let {ai)i<n and {bj)j<n be the diagonal entries of A and B. Observe that 
X{{A + PBS^)SS^{A + PBS^f) = x{{AS + PB){AS + PBf), 

so that 

f[A,B]{x‘^) =Edet(x/-dil(^5 + PS)). 

Let a, TT be the permutations corresponding to S, P, and let ■(— Oi] denote a matrix with 

entry Oj in position ia{i) and zeros elsewhere. Then the nonzero entries of M = AS + PB are 
given by [icr{i) •(— eiOi] + [vr(j)j •(— f,jbj], where and the signs ei,f,j € ±1 are chosen independently. 
In particular, each entry of M has mean zero. 

Thus, if we expand det(x/ — dil(M)) in permutations for any fixed a, vr, the only terms that 
survive are the involutions, which correspond to bipartite matchings in the support of dil(M), 
i.e., pairings of entries from the upper right and lower left blocks of dil(M) (see the appendix 
of [MSS 15a] for the simple calculation). Indexing such /c-matchings of the entries by bijections 
p : S ^ T, where S, T are /s-subsets of the rows and columns of M, we have for fixed every fixed 
(7, vr: 


E,,5det(x/ - dil(M)) = E,,^ n KO)- 

k=0 ^ )-T j^S 

In particular, the coefficient of is: 

5'P-5'^T i^S 

= ^ ^ n(®* l{p(*)=o'(d} + 

p:S—^T i^S 

since the cross terms vanish because of the independent mean zero signs. 

“ ^ ^ ^ ( n ) I n ^p(j)^{^^(*)=d 

g p:S^T RCS \iGR / \ieS\R 

We now calculate the expectation of this with respect to a, vr. The key point is that the events 
{^{p(i)=a{i)}}i€R and the events {l{ 7 rp(i)=i}}iG 5 \_R are independent because a and vr are indepen¬ 
dent and p{R) is disjoint from p{S \ R). Thus, the expectation over cr, tt is given by: 


yz yz i n i{p(i)=cr(i)} 

S^T€{}f)P-S^T Res \ i£R 


Ett fep(j)l{7rp(i)=i} 

ieS\R 


yz y^ y^ (n n ^{pw='7(*)}) i n ^pd ) i n ^{^pw=d 

5:^'7’g([^l) p.S^T ReS \i£R / \ idR / \idS\R J y idS\R 


^In what follows, it is good to think of the matrix M as a directed bipartite graph in which the left nodes 
are indexed by the entries Oi, the right nodes are indexed by the bj, and the entries of M are given by rightward 
edges labeled by aio-(i) and leftward edges labeled by 
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Observe that 


F TT1 _ l-^D- 

li i{p(*)=a(*)} - —— 

ieR 


only depends on the size of R, and similarly 

IEtt l{7rp(j)=i} 

i£S\R 

Thus, the expectation of interest simplifies to; 


dl 


E E E(n»?) ( n 

^ p:S^T R(ZS \iGR / \i^S\R 


d\ 


^ ^ ^ ^^2'^R^i^2\zid—\R\y-{d-\z\y. 


o-ppUd]'! Res p:S^T 
ZCT p(S\R)=Z 
R+Z=k^^ 


d\ d\ 


where we have written Z = p{S \ R) and := a* to ease notation 


i£S 


E E ■ (l-RI)H|2|)r 


5Tg(M) 

\R\ + \Z\=k 


E e E(“i)'*('>: 


r-\-z=k 


r+z=k 


2\Z 

(di)2 

^ ’ \R\=rSeR 

\Z\=z^^Z 

r\z\{d - ry.{d - zy. ^ (d-r\fd-z 

,k — r J \ k — z 




|iJ|=r 

IZI=z 


y. r!zl(d - r)!(d - z)! (d - r)! _ {d - z)'. sy / 

(RU2 (R _ _r-'\Ud- hV(h - zV . 


r-\-z=k 


E 

r+z=k 


2\Z 

(dip (d - ky(k - ry. (d - ky.(k - z)! 

IZI=z 


(d — r)l(d — z)l 

(d-ky.d\ 


er(AA^)ePBB^), 


as desired. 


□ 


3 Real Rootedness of the Asymmetric Additive Convolution 

We will use the theory of stable polynomials to prove Theorem 11.41 For this theorem, we will 
require Hurwitz stable polynomials. We recall that a multivariate polynomial p( 2 ;i,... ,Zm) G 
1R[2;i, ..., Zm] is Hurwitz stable if it is identically zero or if whenever the real part of Zi is positive 
for all i, p(zi,... ,Zm) pO. 

Instead of proving Theorem ll.4l directlv. we prove the following theorem from which it follows. 
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Theorem 3.1. Let 


d d 

p{x) = ^ x'^~'^ai and q{x) = ^ x'^~^bi 
2=0 2=0 

have only nonpositive real roots. Then, 


r{x) x^ ^ 


k=0 


( {d — iV-jd ~ jV- 


2 

aihj 


is also in P (d). 

We will use the following result to prove that a polynomial is in P“. 


Lemma 3.2. Let r{x) he a polynomial such that h{x,y) = r{xy) is Hurwitz stable. Then, 
r E P“. 


Proof. Let ^ be any root of r. If ^ is neither zero or negative, then it has a square root with 
positive real part. Setting both x and y to this square root would contradict the Hurwitz stability 
of h. □ 


We will prove that p l+Pd g is in P“ by constructing a Hurwitz stable polynomial and applying 
Lemma 13.21 To this end, we need a few tools for constructing Hurwitz stable polynomials. 

The first is elementary. 

Claim 3.3. If p{x) E P“, then the polynomial f{x,y) = p{xy) is Hurwitz stable. 

Proof. If both X and y have positive real part, then xy cannot be a nonpositive real, and thus 
cannot be a root of p. □ 

The second is essentially a result of Borcea and Branden. It’s proof is identical to that of 
Proposition 3.4 in [BBDQb] . 

Proposition 3.4 (Polarization). Let 


d d 

V{x-,y) = '^'^Cijx^y^ 
i=0 j=0 

be a Hurwitz stable polynomial. Let for each integer i, let af be the ith elementary symmetric 
polynomial in the variables xi,..., Xd, and let be the corresponding polynomial in yi,..., y^,. 
Then, the polynomial 


P{xi,...,xd,yi,.. 


d d X y 

2 = 0 j=0 \i) \j) 


is Hurwitz stable. 

The polynomial P is called the polarization of p. We remark that P{x,..., x,y,... ,y) = 
P{x,y). 

The last result we need is due to Lieb and Sokal |LS81] (see also |BB09al Theorem 8.4]). 
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Theorem 3.5. Let P{zi,..., Zd) and Q{zi ,..., Zd) he Hurwitz stable polynomials. Let Df denote 
the derivative with respect to Zi. Then, 

Q{Dl,...,Dl)P{zr,...,Zd) 


is Hurwitz stable. 


Proof of Theorem \3.1[ Define f{x, y) = p{xy) and g{x, y) = {xyYq{l/xy). Let F{xi, ..., Xd, yi, • ■ ■ 
be the polarization of f{x, y) in the variables xi, ..., Xd, yi,... ,yd- Let G{xi,..., Xd, yi,..., yd) 
be the analogous polarization of g{x,y). 

Let fjf be the ith. elementary symmetric function in xi,..., Xd, and let Sf be the ith. elemen¬ 
tary symmetric function in Df,..., D^. Define and analogously. 

Then, 

d X V 

F{xi,...,xd,yi,...,yd) = l^ai ——, 

i=0 (j 


and 


Define 


G{Df,...,D-„Dy,... 


^Dy,^±bM. 


d) 


1=0 


cir 


H{xi, ...,Xd,yi,...,yd) = G{Df,..., D^, D\,..., D^)F(xi ,... ,Xd,yi,..., yd)- 
We know from Theorem 13.51 that H is Hurwitz stable. Define 

h{x,y) = H{x,...,x,y,...,y). 

It is immediate that h is Hurwitz stable too. We will prove that h{x,y) = r{xy), which by 
Lemma 13.21 implies that r is in 

It will be convenient to know that 


hi (Tj = 


(‘^+- ifi<i 


otherwise. 


We may now compute 




H{xi,...Xd,yi,...,yd) = 


rr^ rr^ 
fJd-j^d-j 


*=0 IJ 3=0 [jj 

E E 

1=0 j:i<d-j D/ \j) 
i+j<d Uj [jJ 


id\^ (dG 
i+ 3 ^d \j) 


ajbi f d + i — {d — j) 


^d-i-j^d-i-j 


,yd) 
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So, 


(^+Jy 

E -l, \ i ) X y 


i+j<d 

d 


/k\ ^ 


k=0i+j=k (i) (,) 


h{x,y) = X] E 


d\ ik\ \ 2 


m 




Xd-kyd-k_ 


k=0 i-\-j=k 
d 

E E — i—j)i 

k=0i+j=k ^ ^ ’ 

So, r(xy) = /i(x, y) and therefore must have only nonpositive real roots. 


{d - i)\{d - j)\\ d-kd-k 
X y 


□ 


4 Transform bounds 


All of our transform bounds are proved using the following lemma. It allows us to pinch together 
two of the roots of a polynomial without changing the value of the Cauchy transform at a 
particular point. Through judicious use of this lemma, we are able to reduce statements about 
arbitrary polynomials to statements about polynomials with just one root. 

Lemma 4.1 (Pinching). Let a > 0, d > 2, and let p{x) G P(d) have at least two distinct roots. 
Write p{x) = 0^=1 “ -^*) where Ai > A 2 > • • • > and Ai > Xk for some k. Then there exist 

real p and p so that p{x) = p{x) +p{x), where 

p{x) = {x — p)"^ (x — Aj) G P(d) and p{x) = {x — p) (x — Aj) G P(d — 1), 

and 


a. maxroot (C/qp) = maxroot (Uap) = maxroot {UaP), 

b. \i> p> Xk, and 

c. p > Xi- In particular, if d > 3 then p has at least two distinct roots. 
Proof. Let t = maxroot (UaP) and set 


We have chosen p so that 


which implies 


2 

l/(t — Ai) + 1 /{t — Xk) 

2 _ 1 1 
t — p t — Xi t — Xk 

Dp{f) _ Dp{t) _ 1 
p(t) p{t) a ’ 
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and thus maxroot {Uap) = t. Our choice of /r also guarantees that t — is the harmonic average 
of t — Ai and t — Xk- Thus, fi must lie strictly between Ai and Afc, which implies part b. We may 
also conclude that t — fi < (l/2)(2t — (Ai + A^)), which implies that 

^ > (Ai + Afc)/2. (9) 


We have 


p(x) =p(x)-p(x) = ((2^ - (Ai + Afc))® - - AiAfc)) (x-Ai). 


This and inequality ([9]) imply that p{x) G F{d — 1). As Ua is linear, we also have {Uap){t) = 0. 
To finish the proof of part a, we need to show that the maximum root of UaP is less than t. The 
one root of p that is not a root of p is 

(tef p AiAfc 
2p — (Ai + Xk) 


To see that t > p, compute 

(Ai - p){p-Xk) 

2p — Xi — Xk ’ 

which we know is greater than 0 because of ([9]) and the fact that p is between Ai and Xk- This 
completes the proof of part a. 

To prove part c, note that 


{2p — (Ai + Xk)){p — Ai) — — 2Xip + )p — {p — Ai)^ > 0. 


□ 


4.1 Symmetric additive convolution 

Theorem 11.11 tells us that if p{x) G F{d) and q{x) = then p{x) ES^ q{x) = Dp{x). As 

maxroot = {d— l)a, the following lemma may be viewed as a restriction of Theorem II.71 

to the case that q{x) = x‘^~^. 

Lemma 4.2. For a > 0, d > 2, and p G P(d), 

maxroot (UaDp) < maxroot (UaP) — a. (10) 

Proof. Ifp = (x — A)'^, then maxroot (t/ap) = X + da, and maxroot (t/cfiAp) = A + ((i— l)a, giving 
equality in (fTOl) . 

We now prove the rest of the lemma by induction on d, with d = 2 being the base case. For 
a polynomial p in P(d), define 

(j){p) = maxroot (UaP) — a — maxroot {UaDp). 

We will prove that (/)(p) > 0 for all real rooted polynomials p. In particular, for every 72 > 0, we 
prove this for all polynomials p G P(d) whose roots are bounded in abosolute value by R. So, 
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let i? > 0 and let p be a polynomial with roots in [—R, i2] that minimizes (j). As the set [—-R, R]'^ 
is compact, there must be some polynomial at which the minimum is achieved. 

Let /? = maxroot (Uap) — a, (3' = maxroot (UaDp), and assume by way of contradiction that 
/?' > /?. If p has at least two distinct roots, then we may apply Lemma fd.ll to obtain polynomials 
p and p with 

maxroot (UaP) — a = maxroot (Uap) — a = f3. 

Observe that all of the roots of p have absolute value at most R. We will obtain a contradiction 
by showing that maxroot > /3'. 

For d > 3, we use the inductive hypothesis to conclude that 

maxroot {UaDp) < /3, 

and thus {UaDp){l3') >0. Asp = p—p, this implies that {UaD^{j5') < 0, and so maxroot (UaDp) > 
13', a contradiction. For the base case of d = 2, we observe that p has degree 1 and so UaDp 
merely equals the first coefficient of p. As this coefficient is positive, we may again conclude 
that {UaDp){/3') > 0, and thereby obtain a contradiction. □ 

Most of the work in the proof of Theorem 11.71 is devoted to the case in which p and q have 
the same degree. When their degrees are different, we use the following consequence of Theorem 
ll.ll to simplify their symmetric additive convolution. 

Lemma 4.3. For p G P(d) and q G P(/i:) with k < d, 

pS\dq= ^(Dp) lEld-i q- 
d 

Lemma 4.4. For a > 0, p = {x — A)'^ for some real A and q G P(d), 

maxroot {Ua{pSd q)) = maxroot (UaP) + maxroot (Uaq) — cud. 

Proof. From the definition ([T]), we know that p{x) q{x) = q{x — A), and so 

maxroot {Ua{p BEI^ q)) = X + maxroot (Uaq) . 

On the other hand, maxroot {UaP) = A + ad. □ 

Our proof of Theorem 11.71 will be very similar to our proof of Lemma 14.21 

Proof of Theorem [13 We will prove this by induction on d. In the base case of d = 1 we observe 
that 

(x — A) ffli (x — /i) = (x — A — p). 

This gives equality in ([5]), as 

maxroot {Ua{x — A)) = A + a. 

If d > k, then by Lemma 031 we have 

maxroot g)) = maxroot {Ua{DpS\d-i q)), 

< maxroot (UaDp) + maxroot {Uaq) — (d — l)a (by induction) 

< maxroot {UaP) + maxroot {Uaq) — da (by Lemma 02]), 
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as required. 

We now consider the case in which k = d. For a polynomial p in P(d), define 

4>{p) = maxroot (Uap) + maxroot {Uaq) — da — maxroot {Ua{pSd q)) ■ 

As before, we will prove that (j){p) > 0 for all polynomials whose roots are in [—R,R], for 
arbitrary large R. So, let ii > 0 and let p be a polynomial in ¥{d) with roots in [—i?, i?] that 
minimizes (j). If p has just one root, then Lemma 14.41 tells us that ^{p) = 0. Otherwise, let 

/3 = maxroot (f7aP) + maxi'oot — da, and /3'= maxroot O')) j 

and assume by way of contradiction that /3' > /3. 

We may apply Lemma l4.1l to obtain polynomials p and p with 

maxroot (C/op) + maxroot {Uaq) — da = maxroot (UaP) + maxroot {Uaq) — da = (3. 

Observe that all of the roots of p have absolute value at most R. We will obtain a contradiction 
by showing that maxroot {Ua{p'S\d q)) > (3'. 

By the inductive hypothesis, we know that 

maxroot {Ua{pSd q)) < P, 

and thus {Ua{pSd q)){l3') >0. As p = p — p, this implies that {UapSd q){/3') < 0, and so 
maxroot ([/apEBd q) > 13', a contradiction. 

□ 

If one carries out the proof with a little more care, one can show that equality can only be 
achieved when p has just one root or when k = d and q has just one root. 

4.2 Symmetric multiplicative convolution 

We begin by considering the case in which p = {x — A)'^. We then have that 


Thus, 

and 

Lemma 4.5. // A > 0, p{x) 


{w) 


1 + w 
w 


A, 


Sp {w) = A. 

{x — X)^ and q{x) G P^(d), then for all w >0 
{w) = Sp {w) Sg {w). 
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Proof. For p{x) = (x — A)'^, one may use either the definition ([2|) or Theorem 11.51 to compute 


As, ('^^) •^^q{x) ('^)) 


p{x)mdq{x) = X^q{x/\). 


‘^p(a;)[Sldq(a:) (^) ^'^q{x) ('^) • 


□ 


With a little care, one can prove that this is the only case in which equality holds in Theo¬ 
rem [1^ 

The finite multiplicative convolution of polynomials of different degrees may be computed 
by taking the polar derivative with respect to 0 of the polynomial of higher degree. We recall 
that the polar derivative at 0 of a polynomial p of degree d is given by dp — xDp (see [MarGGl 
p. 44]). 

Lemma 4.6. For p{x) G P(d) and q{x) G P(A:), for k < d, 

p{x) [x]^ q{x) = \{xDp{x) - dp{x)) [x]^_i q{x). 
d 

Proof. Follows from Theorem 11.51 bv an elementary computation. □ 


Let R be the operation on degree d polynomials that maps p{x) to x'^p{l/x). The polar 
derivative may be expressed in terms of R by 

dp — xDp = RDRp. 

Claim 4.7. For p = {x — \Y, 

xDp — dp = Xd{x — XY~^. 

For p G P'''(d), {xD — d)p G P'''((i — 1) and 

maxroot (p) > maxroot {xDp — dp). (11) 

Proof. The first part is a simple calculation. Inequality (llip follows from the fact that p G P'’'((i) 
implies that Rp G P^(fi) and the fact that the roots of DRp interlace those of Rp. To see that 
{xD — d)p G P^((i — 1), it remains to check that the lead coefficient is positive. □ 


As we did with symmetric additive convolutions, we relate the A4-transformation to the 
maximum root of a polynomial. We have 

Mp (z) = w maxroot ( ( 1 — Pix)] = 0- 


=> maxroot 
We therefore define the operator Vw by 

Vwp{x) = ( 1 - 


xD 


d{w -|-1) 


d{w -|- 1) 
p{x), 


which gives 

Mp {z) = w maxroot {VwP) = -Z- 

Note that the polar derivative is dVo. 

Our proof of Theorem 11.91 will also employ the following consequence of Lemma 14.11 
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Corollary 4.8. Let w > 0, d > 2, and let p{x) E have at least two distinct roots. Then 

there exist p E and p E F~^{d — 1) so that p{x) = p{x) +p{x), the largest root of p is at 

most the largest root of p, and 

maxroot ( 14 )P) = maxroot (Vy^fP) = maxroot ( 14 )P). 

Proof. To derive this from Lemma l4.ll let t = maxroot {V^p) and set 

tD 

^ d(w + 1) ’ 

The polynomials p and p constructed in Lemma l4.1l now satsify 

maxroot {Uap) = maxroot {UaP) = t = maxroot {Vyjp) = maxroot (V^p) , 
as desired. □ 


Proof of Theorem M.tA As we have already shown that equality holds when one of p or g has 
just one root, we will henceforth assume that both have at least two distinct roots. Under this 
assumption, we will prove that 

w 

maxroot ( 14 >p[x]^ q) < -maxroot (Vyjp) maxroot (Vwq). 

I + w 

We will prove this by induction on the maximum degree of p and q, which we call d. The base 
case of d = 1 follows from Lemma 14.51 

If the degree of p is larger than the degree of q, then we may prove the hypothesis by 


maxroot (14;(plE]d q)) = maxroot {{xDp{x) — dp{x)) q{x)) (by Lemma 


< 


w 


< 


1 + w 
w 

1 + rc 
w 


1 + w 

by applying Claim W7!\ to Vwp{x). 

For polynomials p E P"*", we define 


maxroot {Vy,{xDp{x) — dp{x))) maxroot {Vwq{x)) (by induction) 
maxroot {{xD — d)Vwp{x))) maxroot {Vy,q{x)) (as 14, and Vq commute) 
maxroot {Vwp{x)) maxroot {Vwq{x )), 


4>(p) ———maxroot ( 14 jp) maxroot (Vyjq) — maxroot (VwplEid q) 

1 + w 


We will prove that 4>{p) > 0 for all polynomials p E P"*". As before, we let R be an arbitrary 
positive number and let p be a polynomial minimizing (/>(p) over the space of degree d polynomials 
with roots in [0,i?]. If p has just one real root, then we know from Lemma 14.51 that f){p) = 0. 
If p has two distinct roots, then we apply Corollary 14.81 to obtain polynomials p and p with 
maxroot (14,p) = maxroot (14,p) = maxroot (14jp). We define 
w 

13 =-maxroot iVwp) maxroot iVwq) and (3' = maxroot {Vw(p^d q)) , 

1 + w 


and assume by way of contradiction that (3' > (3. By induction, we know that maxroot (14 (p lElrf q)) < 
(3, and so Vw{p^dQ)iP') > 0. This implies that Vw{p\Eidq){f3') < 0, and thus maxroot {Vw{p [E],^ q)) > 
j3'. As all of the roots of p are less than R, this contradicts the assumption that p minimizes 

m- □ 
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4.3 Asymmetric additive convolution 

The following is a restatement of Theorem 11.81 

Theorem 4.9. Let p{x) and q{x) be in P'’'(4) for d> 1. Then for all a >0, 

maxroot {UaS’ip^^d q)) < maxroot (Ua^’P) + maxroot (C/q,S(?) — 2ad, 
with equality only if p or q equals x^. 

We remark that if q{x) = then p l++l ^ q = P, and 

UaS^q = Uax'^^ = x‘^^~^{x — 2da), 


so 

maxroot {Ua'S>q) = 2ad. 

This is why the theorem holds with equality when q{x) = x'^. 

The following lemma tells us that it suffices to prove Theorem 14.91 in the case that a = 1. 


Lemma 4.10. For a real-rooted polynomial p{x), 

maxroot ([/qp(x)) = —maxroot (t/ip(ax)). 
a 

Proof. Let q{x) = p{ax), so 

Uiq{x) = p{ax) — ap'{ax). 


Let 


w = a-max (p(x)) = maxroot (Uap) 


p{w) — ap'{w) = 0. 


Then, 


{Uiq){w / a) = p{w) — ap\w) = 0. 


□ 


Our proof of Theorem 14.91 will use the following lemma to pinch together roots of p to reduce 
the analysis to a few special cases. 

Corollary 4.11. Let a > 0, d>2, and let p{x) G have at least two distinct roots. Then 

there exist p G and p G P'''(fi — 1) so that p{x) = p{x) +p{x), the largest root of p is at 

most the largest root of p, p has a root larger than 0, and 

maxroot (?7 q,Sp) = maxroot (?7 q§p) = maxroot (?7 q,Sp) (12) 

Proof. Let t = maxroot (C/q,5p), so 

maxroot ((1 — 2atD)p) = y/t. 

Apply Lemma 14.11 with 2atD in the place of a to construct the polynomials p and p. They 
satisfy 

maxroot (Uap) = maxroot {UaP) = Vt, 

which implies (11211 . □ 
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We will build up to the proof of Theorem 14.91 bv first handling three special cases: 


• When p{x) = {x — A)'^ and q{x) = x'^~^. That is, we consider DxD{x — A)'^ (Lemma 14.14p . 

• When p(x) G P^((i) and q{x) = x^~^. That is, we consider DxDp{x) fLemma I4.15p . 

• When p{x) = (x — A)*^ and q{x) = (x — iLemma I4.16p . 

As with the other convolutions, we may compute the asymetric additive convolution of two 
polynomials by first applying an operation to the polynomial of higher degree. In this case it is 
DxD, also known as the “Laguerre Derivitive”. 

Lemma 4.12. Let p € P+((i) and let q G P+(A:) for k < d. Then, 

q = {l/d^){DxDp) ETId-i q. 

Proof. Follows from Theorem II.31 □ 

The following characterization of the Laguerre derivitive also follows from Theorem 11.31 
Claim 4.13. If q{x) = x^~^, then 

p l++lc i q = DxDp. 

Lemma 4.14. For a>0, X>0,d>2 and p{x) = {x — A)'^, 

maxroot [Ua^DxDp) < maxroot {Ua^p) — 2a, 
with equality only if X = 0. 

Proof. By Lemma [4.101 it suffices to consider the case of a = 1. As Sp{x) = — A)*^, 

UiSp{x) = {x^ - 2Xd - X){x^ - X)‘^-\ 

So, the largest root of this polynomial is the largest root of 

rx{x) — 2Ad — A). 

We may also compute 

UiSDxD{x^ - Xy = qx{x){x^ - Xf-\ 

where 

qx{x) dx^ — 2d{d — l)x^ — (d + l)Ax^ + 4(d — l)Ax + A^. 

We now prove that 

maxroot {qx) < maxroot {rx) — 2, 

with equality only if A = 0. 

We first argue that qx{x) is real rooted. This follows from that fact that it is a factor of 
UiE>DxD{x — A)'^. For A > 0 all of the roots of DxD{x — A)'^ are nonnegative, and so applying 
§ to it yields a polynomial with all real roots. 
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We now compute 

maxroot (rx) = d + x/cP + X. 

Define 

fj,x — d 'x/ d‘^ X — 2. 

Elementary algebra gives 

= 4A - 8dx/d^ + X + 8d^ = {2d - 2x/d‘^ + Xf. 

So, qx{nx) > 0, with equality only when A = 0. With just a little more work, we will show that 
fix is an upper bound on the roots of qx for all A. 

For qx to have a root larger than fix, it would have to have two roots larger than fix- When 
A = 0, the polynomial qx has one root at fio and a root at 0 with multiplicity 3. As qx is real 
rooted for all A > 0 and the roots of qx are continuous functions of its coefficients, and thus of 
A, we can conclude that for small A all but one of the roots of qx must be near 0. Thus, for 
sufficiently small A, qx can have at most one root greater than fix, and so it must have none. As 
the largest root of qx and fix are continuous function of A, maxroot (qx) > fix for all sufficiently 
small A. As qxifJ-x) > 0 for all A > 0, we can conclude that maxroot (qx) > fix for all A > 0. □ 

To see that Lemma 14.141 is equivalent Theorem 14.91 in the case of q = , note that for 

q{x) = x^~^, 

UaSq{x) = Uaq{x^) = = x‘^^-^{x - 2{d - l)a). 

The equivalence now follows from Claim 14.131 and the fact that the the largest root of this 
polynomial is 2{d — l)a. 

Lemma 4.15. For a>f),d>2 and p G 

maxroot [Ua^DxDp) < maxroot {Ua^p) — 2a, 

with equality only if p{x) = 

Proof. For every p G P"*", define 

(j){p) = maxroot (UaSp) — maxroot (UaSDxDp) — 2a. 

We will show that (/>(p) >0 for every polynomial p G P"*" of degree at least 2, with equality only 
when p = x^. 

Our proof will be by induction on the degree of p. For any R > 0, let p be the polynomial in 
P^((i) having roots in [0, ii] that minimizes ^{p). Since this set is compact, there is a polynomial 
p at which the minimum is obtained. If all the roots of p are the same, then Lemma [4.141 implies 
that 0(p) > 0, with equality only when all roots of p are zero. If p has two distinct roots, then 
we apply Corollary 14.111 to obtain polynomials p and p. 

Let 


/3 = maxroot (C/oSp) — 2a = maxroot (f7a§p) — 2a = maxroot {UaS'P) — 2a. 
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Now, assume by way of contradiction that 4‘{p) < 0. This means that 

/3' maxroot (Ua^'DxDp) > (5. 

If d > 3, then we may assume by induction that (t>{p) >0 and thus 

Ua^DxDp{l5) > 0. (13) 

If d = 2, then p has degree 1 and so Ua^DxDp equals the lead coefficient of p, which also implies 

m- ^ 

As p = p — p, we may conclude that 

UaSDxDp{l3') = Ua^DxDp{f3') - UaSDxDp{f3') = -Ua§DxDp{/3') < 0, 

and thus 

maxroot {UaS'DxDp) > /3'. 

But, this contracts the assumption that p minimizes (p{p) over polynomials in P'’'(d) with roots 
in [0, R]. □ 

In Section 14.41 we establish the following special case of Theorem 14.91 

Lemma 4.16. For X,p,>0, and d > 1, let p{x) = {x — A)'^ and q{x) = {x — pY- Then for all 
a >0, 

maxroot {UaS{pSSd q)) < maxroot (I7 q,§p) + maxroot {UaEiq) — 2ad. 

We now use Lemma [4.16l to prove Theorem 14.91 through a variation of the pinching argument 
employed in the proof of Lemma 14.141 

Proof of Theorem \4.S\ We will prove this by induction on the maximum degree of p and q, which 
we call d. Our base case of d = 1 is handled by Lemma 14.161 

Assume, without loss of generality, that the degree of p is at least the degree of q. If the 
degree of p is larger than the degree of g, then we may prove the hypothesis by 

maxroot {Ua^ip^Sd q)) = maxroot {Ua^{{DxDp) 1++! ^; i q)) (by Lemma 14.121) 

< maxroot {UaSi{DxDp)) + maxroot [Ua^q) — 2a{d — 1) (by induction) 

< maxroot (t/o-Sp) + maxroot (I7q-§(?) — 2ad (by Lemma l4.14p . 

Lemma l4.141 also tells us that equality is only achieved when p = x^. 

We now consider the case in which both p and q have degree d. For polynomials p and q in 
^■'■(d), define 

4>{P: q) = maxroot (UaSp) + maxroot (UaSq) — 2ad — maxroot (C/Q,S(pim(i q )). 

We will prove that (p{p, q) >0 for all such polynomials, with equality only if one of them equals 
x^. For any R > 0, let p and q be the polynomials with all roots in [0, i?] that minimize (j){p, q). 
Since the set of degree d polynomials with roots in [0, R] is compact, there exist polynomials p 
and q on which the minimum is obtained. If p and q each have only one root, then Lemma 14.161 
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implies ct){p,q) > 0, with equality only if one of them equals x'^. If not, let p have at least two 
distinct roots. We may thus apply Corollary 14.111 to obtain polynomials p and p. 

Let 


/? = maxroot (t/a§p) + maxroot {Ua^q) — 2ad 
= maxroot (t/^Sp) + maxroot {Ua^q) — 2ad 
= maxroot {Ua^p) + maxroot {Ua^q) — 2ad. 

Now, assume by way of contradiction that 4>{p, q) < 0. This means that 

(5' maxroot ([/aS(pimrf q)) > fd. 

For d >2, we may assume by induction that 4>{p, q) > 0 and thus 

(t/„S(pETrfg))(/3) >0. (14) 

As p = p — p, we may conclude that 

([/aS(p^^q))(/d') = ([/aS(p^^q))(/3') - ([/aS(p^^q))(/d') = -([/aS(pES^q))(/d') < 0, 
and thus 

maxroot {UaS'DxDp) > fd'. 

But, this contradicts the assumption that p and q minimize 4>(p, q) over polynomials with roots 
in [0, i?]. 

□ 


4.4 Chebyshev Polynomials 

This section is devoted to the proof of Lemma [4.161 For positive A and p, let 

q^’^{x) = {x-Xf^d{x-p)'^. 

Most of the work will be in the proof of the following lemma. 

Lemma 4.17. 



t 

- (A + //))2 -4A/i 


Proof of Lemma Compute 

U^/dSix - A)-" = U^/dix^ - A)'" = (x^ - Xy-\x^ - X 
So, the maximum root of this polynomial is the maximum root of (x^ 


2q:x). 

A — 2ax), which is 


a 


+ + A. 


Now, define 

t = -x/a^ + A + ■\/ + p. 
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To prove Lemma 14.161 it suffices to show that the maximum root of U^/d^Qd'^ is less than t, 
which follows from Lemma 14.171 bv computing 


(t) < 


2dt 

- (A + /i))2 -4A/X 


d/a. 


□ 


We will bound the Cauchy transform of 
second kind, which we recall are defined by 


by relating it to Chebyshev polynomials of the 


Ld/2J 

Udix) = (2x)''-2^(-1)^ 


k=0 


d-k 

k 


Lemma 4.18. 

Proof. The lemma is an immediate consequence of the following two identities, which we will 


now prove. 

Qd^i^) = 

Ud{x/2 - 1), 

(15) 

and 


/ X (\/A-^)2\ 

1 ^/\fi ^/\fx ] ' 

(16) 


We recall that the Chebyshev polynomials of the second kind may also be defined by the 
recurrence 



Uo{x) = 1 
Ui{x) = 2x 

Ud{x) = 2xUd-i{x) - Ud- 2 {x). 


We now develop an analogous recurrence for the polynomials q^’^{x). We begin by using Theo¬ 
rem Ol to derive 




d 

T.- 

fc =0 


,d—k 


(- 1 ) 


k 

‘E 

i=0 


d — i 
k — i 




By examining the coefficients of * in this expression, we derive 

go’^(x) = 1 

= x - {\ + ^i) 

qYix) = x-{X + fi)qYi^) - ^^Qd-2(^)- 

We may now prove (USD by a simple change of variables. If we substitute y = x/2 — 1 and 
write the recurrence for Ud{y), we obtain exactly the recurrences for q^'^{x). It is particularly 
important to check that this holds for the degree 0 and 1 polynomials in the sequences. 
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We may similarly prove (fT6|) by induction. The base case of degree 0 is trivial and for degree 
1 we have 


1,1 / X ^ ^ \/^)^ 


\VXii 

And, the recurrence for gives 




y/Jql 


-2 = X - (A + /i). 


= ( 


/ X 



y/Xfl 

i "" 


\\/Ah 

y/Xjl 

X 


y/Xfl 

y/Xfl 


- 2 




\A7t \/A^ 

d’M^) - (by induction) 






= (x - (A + /i))5^1^^(x) - 

Proof of Lemma \4-l'A Applying Lemma 14.181 we obtain 

yHnL>^ y) - TT3 „ A,„ (^) 


\,fi 




2d 

2t L)q^'^ 

2d ^ 

2d 

i DUd,,2 


2^/XJld Ud 


((t^ — (A + ^))/2y'^). 


Theorem 14. 191 below tells us that 




(x) < 


1 


d Ud ^ ^/x‘^ — 1 

Substituting this into the above equalities gives 

^ ^ t 1 

(^) ^ 


Theorem 4.19. For t > 1, 


2-\/A^ \/((t^ - (A + ^j))/2y/Xff)‘^ - 1 
t 1 

2-\/A^ ^- (A + ;u))2/4A;U — 1 
t 

- (A + /i))2 -4A// 


(^) < 




□ 


□ 
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We will require a few standard facts about Chebyshev polynomials. Let denote the degree 
d Chebyshev polynomial of the first kind. 

Recall that for x > 1, 


Td{x) 

Ud{x) 


cosh(((i + l)acoshx) 
sinh(((i + l)acoshx) 
sinh acoshx 


where 


cosh(0) = ^ + e , sinh(0) = ^ - e , 


and sinh(acoshx) = \/ — 


We also recall that for x > 1, acoshx > 0. 

We also need the derivative of Ud- 

= 3 - 2^1 + ^)Td+i{x) - xUd{x)) . 

Proof of Theorem\4.19\ Using claim H?20l we compute 




dUd{x) 


1 1 

d x^ — 1 
1 1 
d x^ — 1 




1 


Vx'^ — 1 



Claim 4.20. For x > 1 and d > 0, 

Ud{x) d +1 d +1 

Proof. For 9 = acoshx, we have 


Td+i{x) 

Ud{x) 


coth((d + 1)0) sinh(0), 


(by claim I4.2nh 


□ 


where we recall that 


coth((d + 1)0) 


1 + e-'^id+i)d 
1 _ g-2(d+l)6» • 
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By substituting a 


So, 


29 and t = l/{d + 1) into claim 0211 we learn that 


coth(((i + 1)0) < 



+ 


1 

d + 1 


coth0. 


Td+i{x) 

Ud{x) 




□ 


Claim 4.21. For some a > 0, let 


F{t) 


1 + e-“/* 

1 - e-"/*' 


For all t G (0,1), 

F{t) < tF{l) + (1 - t)F{0) = {l-t)+ tF{l). 

Proof. The above statement with “<” will follow from a demonstration that F is convex on 
(0,1). The inequality is strict because F is not linear. 

The second derivative of F with respect to t is 

(xepl^iflt + 2{a — t)e“/*) 

^4(go/t _ 

We need to show that this is positive for t G (0,1). The only term that is not obviously positive 
is 

2t + 2(a-t)e“/*. 

To show that this is positive, it suffices to show that 

>t-a. 

Dividing both sides by t, this becomes equivalent to 

g-“A > 1 — a/t, 


which is true for all positive a and t. 


□ 


5 Conclusions 

This paper is part of an effort motivated by the method of interlacing polynomials [MSSlSal 
IMSS15b| to better understand how operations that preserve the real rootedness of polynomial 
affect their roots. In |MSS15c| we use the results of this paper to demonstrate the existence of 
bipartite Ramanujan graphs of every number of vertices and every degree. 

The are other interesting connections between the finite free convolutions of polynomials 
introduced here and the standard free probability. For example, Guionnet and Shlyakhtenko 
have proved that in an appropriately taken limit, the inequalities in Theorems 11.71 fL8l and fL9l 
become equalities. We are preparing a paper on this and other relations. 
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